Unusual resonances in nanoplasmonic structures due to nonlocal response 
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We study the nonlocal response of a confined electron gas within the hydrodynamical Drude model. 
We address the question whether plasmonic nanostructures exhibit nonlocal resonances that have 
no counterpart in the local-response Drude model. Avoiding the usual quasi-static approximation, 
we find that such resonances do indeed occur, but only above the plasma frequency. Thus the 
recently found nonlocal resonances at optical frequencies for very small structures, obtained within 
quasi-static approximation, are unphysical. As a specific example we consider nanosized metallic 
cylinders, for which extinction cross sections and field distributions can be calculated analytically. 
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Nanoplasmonics [l|, Q is presently entering an exciting 
era where the metallic structures offer nano-scale features 
that will eventually allow both photons and electrons to 
exhibit their full wave nature. This regime challenges the 
existing theoretical framework resting on a local-response 
picture using bulk-material parameters. In tiny metal- 
lic nanostructures, quantum confinement [3j-|7| and non- 
local response [8l4l8j are believed to change the collec- 
tive plasmonic behaviour with resulting strong optical 
fingerprints and far-reaching consequences for e.g. field- 
enhancement and extinction cross sections. Within non- 
local response, Maxwell's constitutive relation between 
the displacement and the electric fields reads 



D(r,u)=e J dr'e(r,r',u) ■ E{r',u). 



(1) 



The dielectric tensor e(r,r',ui) reduces to e(r,u)5(r — 
r') in the local-response limit. Historically, there has 
been a strong emphasis on nonlocal response in extended 
systems with translational invariance (TI) [10|, where a 
fc-space representation is useful. However, for the present 
problem of metallic nanostructures, TI is broken and a 
real-space description is called for. 

Recent theoretical studies of nanoscale plasmonic 
structures have predicted considerable differences in the 
field distributions and scattering cross sections between 
local and nonlocal response theories, both in numeri- 
cal implementations of a simplified hydrodynamic Drude 
model [l44l8l |. and in corresponding analytical calcu- 
lations [15j. Importantly, new resonances of the free- 
electron plasma were found, also at optical frequencies, 
which have no counterparts in local-response theories. 
Such novel resonances have already gained interest both 
from a fundamental 0] and an applied fl9| perspective. 
At present, the status of these optical nonlocal resonances 
is unclear, since in Ref. [l3| the same nonlocal model was 
used as in Refs. Il4 - |l8l . and yet no corresponding new 
modes were found at visible frequencies. Resolving this 
issue is important for the engineering of ultrasmall plas- 
monic structures with new functionalities 191 21 1. 



In this paper we report that unusual resonances due to 
nonlocal response do exist in nanoplasmonic structures, 
but only above the plasma frequency, not in the visible. 
We illustrate this property of arbitrary plasmonic struc- 
tures by exact calculations for metallic cylinders. We 
also clarify that different implementations of the com- 
mon quasi-static approximation [{J Hi! a re the reason for 
the conflicting results in 14 - 18| and |l3| . Here we refrain 
from making this approximation altogether, and by com- 
parison analyze the validity and implementation of the 
quasi-static approximation in the hydrodynamic model. 

The hydrodynamic Drude model. We express the col- 
lective motion of electrons in an inhomogeneous medium 
in terms of the electron density n(r,t) and the hydro- 
dynamical velocity v(r,t) [ij]. Under the influence of 
macroscopic electromagnetic fields E(r, t) and B(r,t), 
the hydrodynamic model is defined via 10] 



v ■ X7] v = —~/v [E - 

m 



/3 2 

vxB]- — Vn, (2) 
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along with the continuity equation dtn = — V • (nv), ex- 
pressing charge conservation. In the right-hand side of 
Eq. @, the 7-term represents damping, the second term 
is the Lorentz force, while the third term is due to the 
internal kinetic energy of the electron gas, here described 
within the Thomas-Fermi model, with /3 proportional to 
the Fermi velocity Vp. In analogy with hydrodynamics, 
the third term represents a pressure that gives rise to 
a nonlocal dielectric tensor, since energy may be trans- 
ported by other mechanisms than electromagnetic waves. 

We follow the usual approach [ll| to solve Eq. © and 
the continuity equation, by expanding the physical fields 
in a zeroth-order static term, where (e.g., uq is the ho- 
mogeneous static electron density), and a small (by as- 
sumption) first-order dynamic term, thereby linearizing 
the equations. In the frequency domain, we obtain 

/3 2 V [V • J] + uj {u + vy) J = iujuj^eaE, (3a) 

for a homogeneous medium, where J(r) = —enov(r) is 
the current density, and w p is the plasma frequency which 
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also enters the Drude local-response function e(u) = 1 — 
ujp/[u)(uj + i-f)]. We focus on the plasma, leaving out bulk 
interband effects present in real metals that could easily 
be taken into account 1J, l22j , as well as band-bending 
effects at the metal surface. 

The electromagnetic wave equation. The retarded lin- 
earized hydrodynamic model is then fully described by 
Eq. (j3"aj) together with the Maxwell wave equation 



V x V x E 



-E + iuj^ioJ. 



(3b) 



In order to see that these coupled equations §5$ indeed 
describe nonlocal dielectric response, one can in Eq. (|3b[) 
rewrite the current density J as an integral over the 
Green tensor of Eq. (|3aj) and the electric field, whereby 
the nonlocal dielectric tensor of Eq. (JT]) can be identified. 

In a local-response description it is commonplace to 
introduce the quasi-static or curl-free assumption that 
Vx£ = 23J . This well-established approximation lies 
at the heart of most treatments and interpretations of 
electromagnetic wave interactions with sub-wavelength 
structures. Intuitively, one might expect that it can be 
extended to the nonlocal case and indeed several nonlocal 
treatments use this assumption [H, [H, 13- 17|. However, 
as we shall demonstrate, one should proceed with care. 

Three models. In this work we solve the Eqs. di- 
rectly, without further assumptions or approximations. 
We also compare the nonlocal model with two other mod- 
els obtained by further assumptions. The curl-free non- 
local model enforces the condition Vx£ = 0, which with 
Eq. (|3a[) implies that also V x J = in the medium. For 
the differential-operator term in Eq. pap , from now on 
denoted Lj, this has the consequence that V[V-] simpli- 
fies to the Laplace operator V 2 , which gives the model 
used by Ruppin in the context of exciton physics in 1241 . 
and recently in plasmonics by McMahon et al. [14| - |l7| 
and also by ourselves [l8j]. Finally, by assuming Lj = 
in the hydrodynamic treatment (|3a[) . the familiar local 
model is obtained, with J and E related by Ohm's law. 

We assume that the static density of electrons % van- 
ishes outside the metal of volume V, while it is constant 
and equal to the bulk value inside V, thus neglecting 
tunneling effects and inhomogeneous electron distribu- 
tions associated with quantum confinement [H, [|| . As a 
consequence, J = outside V for all three models. 

Boundary conditions. In the local model the current 
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J has the same the spatial dependence as the .E-field. 
Thus, in this case there are no additional boundary con- 
ditions (ABCs) to those already used in Maxwell's equa- 
tions. For the nonlocal-response models on the other 
hand, ABCs are in general needed 13, HI 2^-27|. From 
discussions in the literature it might appear that the 
number of necessary ABCs is a subtle issue, but we em- 
phasize that there should be no ambiguity. The crucial 
point is that the required number of ABCs depends on 
the assumed static electron density profile at the bound- 
For the present problem with the electron 
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aries 

density vanishing identically outside the metal, only one 
ABC is needed in the nonlocal model to obtain unique 
solutions 27J , and it is readily found from the continuity 
equation and Gauss' theorem: h ■ J = on the bound- 
ary, where n is a normal vector to the surface, i.e. the 
normal-component of the current vanishes lfl, 2^, 27 1, 
for all three models. On the other hand, in general the 
tangential current h x J is non-zero. This 'slip' of the 
current is not surprising, since the hydrodynamic equa- 
tion ([2]) describes the plasma as a non- viscous fluid. 

Likewise, in several implementations of the quasi-static 
approximation, no further ABCs are needed to uniq uely 



determine the electric field and current density [11 
In contrast, in the curl- free nonlocal model of Refs 
24j one more ABC is needed 
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Assumed is that 
the tangential components of J vanish at the bound- 
ary (n x J = 0), so that both normal and tangential 
components of the current field vanish on the boundary. 
In the different context of exciton physics [HI these are 
often referred to as Pekar's additional boundary condi- 
tions. There, the vanishing of the tangential boundary 
currents is motivated by the physical assumption that ex- 
citon wave functions vanish on the boundary [24l |28| . In- 
stead, in the hydrodynamical theory of metals, the ABC 
n x J = seems more ad hoc: not a direct consequence 
of the quasi-static approximation, and not correct if that 
approximation is not made. The different boundary con- 
ditions are summarized in Table HI 

Extinction cross section of metallic nanowires. To il- 
lustrate the surprisingly different physical consequences 
of the three models, we consider light scattering by a 
nanowire. Rather than solving Eqs. @ numerically for a 
general cross-sectional geometry, we here limit our anal- 
ysis to cylindrical wires whereby significant analytical 
progress is possible. We use an extended Mie theory, 
developed by Ruppin 0, Hil , to calculate the extinction 
cross section <j ex t of an infinitely long spatially dispersive 
cylindrical metal nanowire in vacuum. Outside the wire 
there are incoming and scattered fields (both divergence- 
free), whereas inside the wire both divergence-free and 
curl-free modes can be excited, the latter type only in 
case of nonlocal response. The cross section is [30( 



TABLE I: Summary of the three different response models. 
V is the volume of the nanostructure, and dV its boundary. 
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where a is the radius, fco = uj/c is the vacuum wave vec- 
tor, and a n is a cylindrical Bessel-function expansion co- 
efficient for the scattered fields. We consider a normally 
incident plane wave with the electric-field polarization 
perpendicular to the cylinder axis (TM). The expression 
for the coefficients a n depends on the particular response 
model and the associated ABCs. For the curl-free non- 
local model, the a n are known [HI- For the full hydro- 
dynamic model we follow the approach of Ref. 0, where 
the ABC of Ref. |26| is employed. This ABC is for metals 
in free space equivalent to n ■ J = 0. We obtain 

a _ [d n + J' n (K t aj\ J n (k a) - y/eJ n {K t a) J»(fc a) 
[d n + J' n {K t a)] H n (k a) - ^/eJ n (n t a)H^(k a a) ' 

(5) 

where J n and H n are Bessel and Hankel functions of the 
first kind and n 2 = e(cj)t^. The d n coefficients are 



n 2 J n {nia) Jn(n t a) 
Kia Jn(ma) K t a 



[s{w) 1] 



(6) 



where nf — (lo 2 + iuj^/ — u) 2 )/ (3 2 . In the limit /3 — > 0, the 
d n vanish and the a n of Eq. © reduce to the local Drude 
scattering coefficients [3(j, which confirms that nonlocal 
response in our model requires moving charges. 

Are there nonlocal resonances? Figure [T] depicts the 
extinction cross section of Eq. (j4|) for two cylinder radii, 
comparing the nonlocal models with the local Drude 
model. The main surface-plasmon resonance peak at 
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FIG. 1: (Color online) Extinction cross sections <j cxt as a 
function of frequency for TM-polarized light normally incident 
on a metallic cylinder in vacuum. Parameters for Au as in 



tjypi is blueshifted as compared to the local model, Ref. B hcj p = 8.812 eV, fry = 0.0752 eV, and vf = 1.39* 



and more so for smaller radii. Similar blueshifts have 
been reported for other geometries 



10 m/s. Inset: frequency shift of the maximum <7 oxt (cj) for 
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121 and in the curl- nonlocal against local response, as a function of radius. 



free nonlocal model 

Figure 1 shows the unusual resonances mentioned in 
the title of this paper: additional peaks do appear in the 
nonlocal theory but only for frequencies above the plasma 
frequency w p (huj p = 8.9 eV for Ag and Au; 1.5 to 3eV 
is visible). These peaks (such as P2 in Fig. [T]) are due 
to the excitation of confined longitudinal modes, which 
are bulk-plasmon states with discrete energies above hcu p 
due to confinement in the cylinder 13j. These peaks are 
analogous to discrete absorption lines above the band 
gap in quantum-confined semiconductor structures. In- 
terestingly, contrary to the common belief that light does 
not scatter off bulk plasmons, which is correct in the lo- 
cal theory (i.e. no peak around lu p in Fig.[lJ, here in the 
nonlocal model we do find such a coupling to longitudinal 
modes. The new resonances could therefore be observed 
with electron loss spectroscopy but also with extreme UV 
light. The curl-free model also exhibits these resonances. 

The striking difference between the two nonlocal- 
response models is that the curl-free nonlocal model 
shows additional stronger resonances, both above and 
below the plasma frequency, such as PI in Fig. [TJ in par- 
ticular also at optical frequencies. These peaks do not 
show up in the full hydrodynamical model, and thus orig- 
inate from a mathematical approximation rather than 



a physical mechanism. It would however be premature 
to conclude that the quasi-static approximation breaks 
down, because in Ref. [l_3| the modes of cylinders in the 
hydrodynamical Drude model were found after making 
the quasi-static approximation, and the only novel modes 
found were the confined bulk plasmon modes above uj p . 
Fig. [l] also illustrates that for increasing radii, er cx t in the 
two nonlocal models converges towards the local-response 
value. This convergence is slower for the curl-free model. 

In Figure HJa) we depict the scaled displacement- 
field distributions for the three models at the frequency 
marked PI in Fig. 1, where only the curl-free nonlocal 
model has a (spurious) resonance. Correspondingly, in 
Fig. EJa) we find a standing-wave pattern only in that 
model. Its appearance in the displacement field illus- 
trates that the spurious resonance is a transverse reso- 
nance, i.e. occurring in the divergence- free components 
of E and J. Fig. [2jb) on the other hand shows the nor- 
malized eZeciric-field intensity for a true resonant mode 
at the frequency P2 of Fig. [1] Only the two nonlocal 
models give rise to resonant electric-field patterns. These 
confined bulk plasmon modes are longitudinal and would 
not produce standing waves in the displacement field. 
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mation that agree with the quasi-static modes of Ref. Il3, 
so we do not claim a general breakdown of the approx- 
imation itself. Even though there are no nonlocal reso- 
nances in the visible, plasmonic field enhancements are 
affected by nonlocal response. For arbitrary geometries, 
numerical methods must be used to quantitatively assess 
their importance. Self-consistent versions of the versa- 
tile time-domain [l4T - [l7| and frequency-domain [l8| im- 
plementations of the hydrodynamical model can do just 
that. 

This work was financially supported by Danish Re- 
search Council for Technology and Production Sciences 
(Grant No. 274-07-0080), and by the FiDiPro program 
of the Finnish Academy. 



FIG. 2: (Color online) Field distributions in the three different 
models, for TM-polarized light normally incident on a cylin- 
der of radius a — 2nm. (a) Normalized displacement field 
\D\ 2 /\D in \ 2 at the frequency u> = 0.6503w p (PI in Fig. [TJ. 
Din = eo-Ein and Bi n is the incident electric field, (b) Anal- 
ogous plots of \E\ 2 /\E- lYl \ 2 for u = 1.1963oj p (P2 in Fig. [[J. 



Origin of spurious resonances. By eliminating the elec- 
tric field from Eqs. (J3J), it follows that the exact hydrody- 
namic current satisfies the pair of third-order equations 



(/3 2 V 2 +lj 2 + iuj-f — ojp) V • J 



fc 2 V 2 



2 e(w)) V x J = 0, 



(7a) 
(7b) 



which reduce to the more symmetric Boardman equa- 
tions 31( in the absence of damping. For arbitrary ge- 
ometry, Eq. (jTa|) has damped solutions of V- J for v < ui p 
and finite- width resonances for uj > w p , as seen in Fig. [T] 
Both solutions can be consistent with the quasi-static ap- 
proximation V x J = that trivially solves Eq. (|7b|) . On 
the other hand, we find that the spurious resonances have 
resonant divergence- free components of E and J. How- 
ever, these cannot at the same time be curl-free. Thus 
the curl-free nonlocal model has resonant solutions with 
nonvanishing curl, which is logically inconsistent. But 
how could this arise? Once the V x J = assumption 
has been invoked to simplify the differential operator into 
Lj = (3 2 V 2 , the resulting Laplacian equation analogous 



to pap carries no information that the resulting solution 
should also be curl- free. Thus, the solutions found for 
this equation are not necessarily self-consistent. 

Conclusions. We have shown that plasmonic nanos- 
tructures exhibit novel resonances due to nonlocal re- 
sponse in the hydrodynamic Drude model, but only 
above the plasma frequency. The recently reported non- 
local resonances in the visible 14 -[H| agree with older 
work |24j , but are a surprisingly pronounced consequence 
of an implementation of the quasi-static approximation 
that is not self-consistent. For nanowires, we find extinc- 
tion resonances without making the quasi-static approxi- 
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